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ABSTRACT 

Novae are cataclysmic variables driven by accretion of H-rich material onto a 
white-dwarf (WD) star from its low-mass main-sequence binary companion. New 
time-domain observational capabilities, such as the Palomar Transient Factory 
and Pan-STARRS, have revealed a diversity of their behaviour that should be 
theoretically addressed. Nova outbursts depend sensitively on nuclear physics 
data, and more readily available nova simulations are needed in order to effec- 
tively prioritize experimental effort in nuclear astrophysics. In this paper we use 
the MESA stellar evolution code to construct multicycle nova evolution sequences 
with CO WD cores. We explore a range of WD masses and accretion rates as well 
as the effect of different cooling times before the onset of accretion. In addition, 
we study the dependence on the elemental abundance distribution of accreted 
material and convective boundary mixing at the core-envelope interface. Models 
with such convective boundary mixing display an enrichment of the accreted en- 
velope with C and O from the underlying white dwarf that is commensurate with 
observations. We compare our results with the previous work and investigate a 
new scenario for novae with the 3 He-triggered convection. 

Subject headings: methods: numerical — stars: novae — stars: abundances — 
stars: evolution — stars: interiors 
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Introduction 



Classical novae are the result of thermonuclear explosions of H-rich material occurring on 
the surfaces of white dwarfs (WDs) . Such events typically occur in a close binary system con- 
taining a cold WD pr imary component and a low-mass main-se quence (MS) star, the latter 
filling its Roche lobe (IJose fc Hernanzll2007at iGehrz et al.lll998f ). The solar-composition ma- 
terial from the envelope of the secondary component streams through the inner Lagrangian 
point to form an accretion disk and eventually, after having lost its orbital angular mo- 
mentum, joins the WD. For sufficiently low accretion rates (M rj 10~ n - 10~ 9 M©/yr), the 
accreted material accumulates in a thin layer atop the WD until its base temperature, that 
rises because of gravitational compression, reaches a value at which H begins to burn, ini- 
tially in the pp-chain reactions and then in the CNO cycle. As a result, a thermonuclear 
runaway (TNR) ensues causing rapid increases in both the temperature and energy output. 
Peak temperatures during classical nova outbursts can be as high as T max = 2 — 4 x 10 8 K, 
approaching the virial temperature. Under such conditions, the proton-capture reactions of 
the CNO cycle become so fast that they build up large amounts of /3 + -unstable isotopes of N, 
O, and F. T heir temperature-independent deca y rates then limit the energy generation in the 
CNO-cycle ( jStarrfield. Truran. fc Sparks! Il978l ). For neon novae, the activation of the NeNa 
and MgAl cycles driven by injection of Ne-seed nuclei from the ONe-rich substrate can lead 
to nuclear processing of even higher mass isotopes. After a time period of 10 2 — 10 3 seconds, 
the electron-degenerate conditions are lifted as the hot convective envelope expands. Mass 
loss ensues either from a wind or Roche lobe overflow triggered by the radius expansion of 
the burning WD. 



Despite much progress over the past decades (e.g. Uose fc Hernanzll2007al ). some key as- 
pects of classical novae are still poorly understood. An important one is the mixing between 
the accreted envelope and WD, which is required to explain the enrichment observed in most 
nova ejecta in heavy elements, such as C, N, O, and Ne. These can reach combined total 
ejected mass fractions of 30-40% (e.g., IGehrz et al.l Il998l ). The TNR peak temperatures 
and durations do not allow for fuelling production, so the only possible explanation of this 
enrichment is that it originates from the underlying CO or ONe WDs. The proposed mixing 
models can be divided into two groups. Either it is assumed that the en yelope and WD get 



mixed in a thin layer close to their interface before the TNR ensues (e.g., iPrialnik fc Kovetz 
19841 ; lMacDonaldlll983t lAlexakis et al.ll2004j ). or mixing is assumed to be the result of hydro- 
dyna mic boundary mixing at t he bottom of the convection zone triggered by the TNR itself 



2010, 2011a 



Glasner. fc Livnd Il995l ; Idasner. Livne. fc Truranl 119971 . 120051 120071 ; ICasanova et al. 



The mechanism of mass loss by novae is also not fully understood. Apparently, a su- 
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personic outflow, a wind driven by the super-Eddington luminosity ( IKato fc Hachisulll994l ). 
and an expansion of nova eject a as a result of its common-envelope interaction with the 
companion star (e.g., iLivio et al.lll990l ) all contribute to the mass ejection. Finally, there 
are a few key nuclear reactions relevant for nova nucleosynt hesis whose rates still re main 



uncertain, e.g. 18 F(p,a), 25 Al(p,7), and 30 P(p,7), according to I Jose fc Hernana (j2007al ). 



Nova outbursts are recurrent events, where depending on the parameters of a particular 
system, 10 3 -10 4 H-shell flashes would occur in a single binary system. However, most 
computer simulations of nova follow only an i ndividual outburst . An exception are the model s 
by the Tel Aviv group (jPrialnik &: Kovetzl Il995l ; lYaron et al.l 120051 ; lEpelstain et al.l 120071 ). 
who followed more than 1000 flashes. They adopted mixing of the first type mentioned above, 
i.e. diffusive interface mixing during quiet accretion phases between consecutive outbursts. 

In this paper we address and investigate some of these issues. We describe the simulation 
assumptions in Section [2j As a start, we have constructed a grid of multicycle simulations 
of CO novae with the masses 0.65 M , 0.85 M , 1.0 M , 1.15 M , and 1.2 M that accrete 
the solar-composition and CO-enriched material at rates M = 10~ n , 10~ 10 , and lO _9 M /yr 
(Section [3]). The exact value of the WD mass dividing CO and ONe novae is not well 
settled. Current es timates suggest a va lue of about 1.1 M® when the effect of binarity is 
taken into a ccount (IGil- Pons et al.l 120031) . As a verification exercise, we confirm the result 
obtained by Glasner fc Truranl ( 20091 ) that, for massive and slowly accreting CO WDs, the 
peak temperature achieved during the TNR becomes a much steeper function of WD's central 
temperature, when the latter is lower than 10 7 K (Section Hj). We then investigate, for the 
1.2 Mq case, the effect of convective boundary mixing (CBM) during the TNR (Section 15. ip . 
For these simulations we invoke an exponentially decaying mixing efficiency model that has 
already been studied previously in related stellar evolu tion phases, such as the He-shell flash 
convection zones in AGB stars (e.g. iHerwig et al.lll999l ). Finally, we report on a new scenario 
for novae with the 3 He-triggered convection (Section 15. 2p . and give conclusions in Section EJ 



2. Simulation Assumptions 



The calculations are per formed with the ste llar evolution code Modules for Experiments 

[rev. 3611] . From the a vailable options we 



in Stellar Evolution (MESA. lPaxton et al.ll2011 



use the 2005 update of the OPAL EOS tables (IRogers fc Nayfonovll2002l ) supplemented for 
lower temperatures and densities by the SCVH EOS that includes partial dissociatio n and 
ionization caused by pressure and temperature (jSaumon. Chabrier. fc van Hornlll995l ). Ad- 



MESA web page with detailed instructions: http://mesa.sourceforge.net. 
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ditionally, the HELM flTimmes fc Swestvl I2OO0I ) and PC flPotekhin k Chabrierll2010[ ) EOSs 
are used to cover the regions where the first two EOSs are not applicable. In particular, 
the HELM EOS takes into account electron-positron pairs at high temperature, while the 
PC EOS incorporates crystallization at low temperature. Both assume complete ionization. 
There are smooth transitions between the four EOS tables. When creating WD models and 
following their subsequent nova evolution, different parts of our computed stellar models are 
entering the p-T domains co vered by all of the above EOSs. We u se the OPAL opacitie s 
(jlglesias fc Rogerslll993l 119961 ) supplement ed by the low T opac ities of Ferguson et al. 



and b y the electron conduction opacities of ICassisi et al.l ( 120071 ) (for details, see 
201lh . 



Paxton et al. 



As in any other phase of stellar evolution, the adopted nuclear network has to be as large 
as necessary in order to account for the energy generation, yet as small as possible in order 
to make computations not too expensive. In MESA, we use the option to solve the nuclear 
reaction network, the structure equation and mixing operators simultaneously. This leads to 
a more stable numerical behaviour, however it makes adding species to the nuclear network 
relatively more expensive compared to using a nuclear reaction operator split option. 

We have started our CO nova simulations with 77 isotopes from H to 40 Ca coupled by 
442 reactions and then gradually reduced these numbers checking that this does not lead to 
noticeable changes in the time variation of the peak temperature. As a result, an acceptable 
compromise has been found empirically, and the following 33 isotopes were selected: 1 H, 
3 He, 4 He, 7 Li, 7 Be, 8 B, n B, 12 C, 13 C, 13 N, 14 N, 15 N, 14 0, 15 0, 16 0, 17 0, 18 0, 17 F, 18 F, 19 F, 
18 Ne, 19 Ne, 20 Ne, 21 Ne, 22 Ne, 21 Na, 22 Na, 23 Na, 22 Mg, 23 Mg, 24 Mg, 25 Mg, and 26 Mg. These 
are coupled by 65 reactions, including those of the pp chains (the pep reaction, whose im- 
portance was emphasized by IStarrfield et al.l (120091 ). has also been added), CNO and NeNa 
cycles. The underlying CO WD models have been prepared using the same isotopes, while 
the reaction list was ext ended to take into account He a nd C burning . By d efault, MESA 
uses reaction rates from ICaughlan. fc Fowlerl (119881 ) and lAngulo et al.l (119991 ) . with prefer- 
ence given to the second source (NACRE). It includes updates to the NA CRE rates for 



14 



N(p,7) 15 Q 



(IGorres et al 



Imbriani et al.ll2005l ). the triple-a reaction (IFynbo et al.ll2005l ). 14 N(a,7) 18 F 
|2000|), and~^C(a,7) 16 faunz et al.l l2002h . Although the main nuclear path 



for a classical nova is driven by p-capture reactions and /3-decays, the a-reactions are im- 
portant for establishing the chemical composition of its underlying WD. As a test, we have 
also tried the MESA s econd option for ch oosing the reactions rates that gives preference to 
the JINA REACLIB (jCyburt et al.l 120101) . This datab ase has many nuclear reaction rates 
in common with those considered by llliadis et al.l ( 120 101 ) . We have not found any significant 
differences in our nova T max (t) profiles between the two options, therefore we decided to stick 
to the default one. 
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The options that specify the physics and numerics assumptions of a MESA simulation 
are set in an inlist file. We have started from inlists of two MESA test suite cases relevant 
for our problem: 

• make_co_wd combines some "stellar engineering" tricks into a procedure that creates 
CO WD models from a range of initial masses (see below), 

• wd2 demonstrates the use of parameters that control accretion, as well as mass ejection 
options relevant for nova calculations. 

Test suite case make_co_wd uses limits on the opacity during the AGB evolution to 
achieve rapid envelope removal without having to compute the details of many thermal pulses 
which eventually culminate in the so-called superwind phase. The simulation of this final 
tip- AGB phase of evolution is not straightforward and beyond the scope of this investigation . 



The procedure followed here is different from that described by lWagenhuber fc Weissl ( 119941 ). 
but has in the end the same effect. 

First, a sufficiently massive initial model has to be chosen, e.g. a 6M pre-MS star for 
the 0.85 M Q WD. Its evolution is computed until the mass of its He-exhausted core reaches 
a value close to the final WD's mass. This phase is shown with the solid blue curves in 
Fig.CD After that, the maximum opacity is reduced to a small value (the red curves), the 
total mass of the star is relaxed to Mwd (the solid black and dashed green curves), the 
maximum opacity is restored (the small red circles), and finally the WD model is given 
time to relax (the dot-dashed cyan curves). As a result, the stellar model arrives at the 
WD cooling track, from which a WD model with a required central temperature, Twd, (a 
lower Twd corresponds to a longer WD's cooling time) can be selected (the dotted magenta 
curves) . 

In reality, a CO or an ONe WD recently formed in a binary system should be surrounded 
by a buffer zone of unburnt material (He-rich or CO-rich, resp ectively), and qu ite a large 



number of nova outbursts have to occur before it will be removed (IJose et al.ll2003l ). Our WD 
making procedure does have a step on which the WD models possess such buffer zones. To 
avoid a discussion of effects to which the presence of He-rich buffer zones can lead, we remove 
them artificially and use naked CO WDs as the initial models for our nova simulations in 
this paper. 

A set of the more massive CO WD models was obtained differently, by letting the 
1.0 M Q CO WD accrete its core composition material until the total masses of 1.15 M & and 
1.2 M e were accumulated, after which the stars were cooled off to generate the massive CO 
WD models for a range of initial core temperature (Tabled]). This special procedure has 
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been employed because it is impossible to obtain a CO WD with Mwd ^ 1.0 M following 
the evolu tion of a massive AG B star, unless one artificially turns off carbon burning. For 
example, iGil-Pons et al.l (120031 ) have found that for M Z ams increasing from 9.3 M to 11 M 
the final WD core mass for the binary star evolution is changing from 1.07 M to 1.22 M 
but, because of carbon burning in the core, it actually contains an ONe WD surrounded by 
a CO buffer zone, the relative mass of the latter decreasing from 7% to 0.8% for the given 
mass interval. Therefore, our massive CO WD models should be considered as substitutes 
for such ONe WDs with CO surface buffer zones. Besides, they are used for comparison with 
other nova simulations that involved CO WDs of similar masses. 



3. MESA Models of Multicycle CO Nova Outbursts 

For either CO or ONe WD the main properties of its nova outburst (such as total 
accreted and ejected masses, peak temperature, maximum luminosity, envelope expansion 
velocity, and chemical composition of the ejecta) depend mainly on the following four pa- 
rameters: the WD mass Mwd, its central temperature Twd (o r luminosity) , the accre- 



tion rate M, and the metallicity of the ac creted material (e.g., iPrialnik fe Kovetzl 11995 



Townsley fc Bildstenl 12004] ; I Jose et al.l 120071 ). In this paper, we consider only CO WDs in 
binary systems with solar-metallicity companions, and in this section we ignore any mixing 
between the WD and its accreted envelope. 

Fig.[2] shows a snapshot of the evolution of the hottest of our 1.15 M CO WD models 
(Twd ~ 3 x 10 7 K; for the correspondence between the WD's initial central temperature 
Twd and luminosity L^d > see Table [T]) accreting the solar-composition material with the 
rate M = 10~ 10 M /yr. The total mass M = Mwd + M acc increases during an accretion 
phase and decreases during a mass-loss event. The abrupt changes in the abundances at 
xwd ~ —4.6 and xqe ~ —8.8 mark the WD's surface and upper boundary of the convective 
envelope, respectively. During an accretion phase, xwd moves to the left, which signifies 
that the mass (of the envelope) to the right of iwd increases; xwd shifts to the right during 
a mass-loss event J§ 

Fig. [2] corresponds to the moment, when the temperature at the interface between the 
WD and accreted envelope has reached its maximum value (a sharp peak on the solid red 
curve in the upper-right panel), and most of the initially abundant 12 C, 13 C, 14 N, and 16 
nuclei have been transformed into the /3 + -unstable p-capture product isotopes 13 N, 14 0, 15 0, 



2 A movie demonstrating the multicycle evolution of the 1.15 M© CO nova model is available at 
http : //astro . triumf . ca/nova-movies. 
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and 17 F (the dot-dashed blue, dashed black, dotted red and solid magenta curves in the 
middle-right panel) in the convective envelope (a region where V ra d > V > V a d in the 
lower-right panel). The /3-decay times of these isotopes limit the energy production, unless 
additional C and O is mixed from the layers below the convection zone. 

For the multicycle simulations mass loss is triggered when the models reach super- 
Eddington luminosities according to the following prescription: 

(L — -^Edd) 



M = -2r] Edd 



1 



where v esc = >J2GM/R, and LEdd = (AnGcM) j 'k. Here, M, R, and L > Lsdd are the mass, 
radius, and luminosity of the star, while k is the Rosseland mean opacity at the surface. The 
scaling factor has been set to ??Edd = 1- This prescription simply assumes that the excess 
of nova luminosity over the Eddington one determines the rate of change of the mass-loss 
kinetic energy. 

The main results of our 1.15 M & CO nova simulations for the mass accretion rate M = 
lO~ lo M /yr are presented in Fig.[3J Cases for accretion of solar and 30% CO enriched 
material, as well as for different WD central temperatures are shown. The case with CO- 
enriched accretion material is not realistic because in the majority of actual CO novae the 
donor star is a Pop I MS dwarf providing material close to solar to be accreted by the WD. 
Therefore, the CO enrichment of nova ejecta comes from the underlying CO WDs and occurs 
either before or during the TNR. Nevertheless, the case with CO enhanced accretor material 
is frequently considered in the lit erature as an artificial way to mimic the effe ct of mixing at 
the core-envelope interface (e.g.. |Jose fc Hernanzlll998l ; IStarrfield et al.lll998l ). 



The effect of the initial WD luminosity (or central temperature) on the strength of 



the nova explosion has been analyzed before (e.g., IStarrfield et al.lll998t lYaron et al. 



2005 



Jose Hernanzll2007bl ). As expected, the strongest outburst occurs in the case of the coolest 
WD (the dotted red curves). Its lower temperature allows the WD to accumulate a slightly 
more massive H-rich envelope on a longer timescale (the lower-right panel) before the TNR 
is triggered. This leads to a higher peak temperature, T max w 2.29 x 10 8 K, and to a longer 
time of envelope's removal by the mass-loss that extends the nova evolution track towards 
lower effective temperatures and larger radii. The corresponding track (dotted red) reaches 
a maximum radius R ~ 3R Q . Evidently, to produce a self-consistent model, we should not 
allow the nova to expand far beyond its Roche lobe radius -Rrl- For example, if the 1.15 M & 
WD has a 0.6 M & MS companion then, for the latter to fill its Roche lobe and therefore be 
able to transfer its mass onto the WD, the binary rotation period has to be nearly 5 hours 
(for a circular orbit with a semi-major axis a = 1.8 Rq). In this case, the WD itself will have 
a Roche lobe with R^ w 0.8 Rq. The MESA stellar evolution code has an option to limit 
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the growth of a star beyond its Roche lobe radius by exponentially increasing the mass-loss 
rate when R > -Rrl, however we did not implement this option in these simulations. 

The CO-enriched material ignites much earlier than the solar-composition case because 
it contains much larger mass fractions of CO isotopes that serve as catalysts for H burning in 
the CNO cycle. As a result of the lower accreted mass, the TNR peak temperature reaches 
only 1.59 x 10 8 K in this case (the dashed blue curves in the lower panels of Fig.|3]). Tabled] 
summarizes the accreted masses, maximum H-burning luminosities, and peak temperatures 
as functions of Mwd> ^wd, -^wd, and M obtained in the nova simulations. The number of 
grid zones in the envelopes of the MESA nova models varies between 500 and 1000, depending 
on the complexity and evolutionary phase of the model. A comparable number of grid zones 
is allocated to the underlying WD. 



4. Comparison with Other Nova Simulations 



Ami Glasner provided us with the parameters of his ID model of a nova outburst 
occurring on a 1.147 M© CO WD. It has a central temperature Twd = 2.4 x 10 7 K, accretes 
solar-com position material at a rate M = 1O~ 1O M / yr, and no interface mixing has been 
adopted (IGlasner. Livne. fc Truranll201lf ). The temperature and luminosity comparison of 
that model with our 1.15 M CO nova simulation with Twd ~ 2.5 x 10 7 K shows very good 
agreement, in particular with respect to the amplitudes (Fig.HJ). The relative differences 
between the accreted masses and peak temperatures for the two models are only 16% and 
6%, respectively. 

As a second test , we c heck if our nova model agrees with findings reported recently 
by IGlasner fc Truranl (120091 ) that the peak temperature achieved during the TNR becomes 
a much steeper function of Twd when th e latter is lower than 10 7 K . Their study was 
motivated by the result obtained earlier by lTownsley fc Bildstenl ( 120041 ). according to which 
such cold WDs should be associated with nova outbursts occurring in binaries with M < 
10~ 10 M Q / yr. To carry out this test, we let the 1.2 M CO WD model cool down to a central 
temperature Twd = 3.3 x 10 6 K, let it accrete the solar-composition material with the rate 
M = 10- n M o /yr, and followed the ensuing nova outburst. This simulation has been 
complemented with the ones done for Twd — 7, 15, 20 and 30 x 10 6 K and the same M. The 
resulting relation between T max and T W d is very similar to those plotted by Glasner & Truran 
in their Fig. 1, and our model therefore confirm their findings. Quantitatively, for Mwd = 
1.2 Mq, T W d = 3.3 x 10 7 K, and M = I0~ u M / yr our model accretes M acc = 1.42 x 10~ 4 M 
envelope mass before TNR ignition and its peak temperature is T max = 3.28 x 10 8 K, while 
Glasner & Truran find M acc = 1.30 x 10~ 4 M Q and T max = 3.68 x 10 s K for their slightly 
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more massive WD (Mwd = 1.25 M ) with Twd = 4 x 10 7 K. Our data are also i n a go od 
agreement with the estimates of M acc and T WD presented by lTownsley &: Bildstenl (12004 ) in 
their Fig. 8. 

Finally, Fig. [5] shows solar-scaled mass-averaged abundances in the expanding envelope 
of a last model from our simulations of a CO nova with Mwd = 115 M , Twd = 15 x 10 6 K, 
and M = 2 x 10~ 10 M Q /yr. These simulations have used an extended nuclear network that 
included 48 isotopes from H to 40 Ca coupled by 120 reactions. The accreted material was 
assumed to be a mixture of 50% solar and 50% W D's core compositions. T his nova model has 
parameters similar to those of the model C05 of I Jose fc Hernanzl (I1998I ). A comparison of 
our fin al abundances from Fig.[5]with those for the model C05 presented by I Jose Hernanz 
(119981 ) in their Fig. 1 shows a very good qualitative agreement. 



5. Effects of Convective Boundary Mixing 
5.1. A Standard Mixing Model 



The nova models presented in Section [3] do not reproduce the ob served enrichment of 



nova ejecta in C, N, O, and other heavy elements (e.g. lGehrz et al.lll998l ) because they do not 
include the interface mixing between the accreted H-rich envelope and its underlying CO WD. 
Recent two- and three dimensional nuclear-hydro dynamic simulations of a nova outburst 
have shown that a possible mechanism of this mixing are the hydrodynamic instabilities 
and shear-flow turbulence induced by steep horizontal ve locity gradients at the bottom of 
the convection zone of the TNR (ICasanova et al.l l2011bl ). These hydrodynamic processes 
associated with the convective boundary lead to convective boundary mixing (CBM) at the 
base of the accreted envelope into the outer layers of the WD. As a result, CO-rich (or 
ONe-rich) material is dredged-up during the TNR. 

Our nova simulations have been performed with the one-dimensional stellar evolution 
code MESA. For one-dimensional CBM calculations, MESA provides a simple model that 
treats the time-dependent mixing as a diffusion process, and that approximates the rate 
of mixing by an exponentially decreasing function of a distance from the formal convective 
boundary, 



D ov = D exp 



2|r — tq\ 
fH P 



(2) 



where H P is the pressure scale height, and D is a diffusion coefficient, calculated using a 
mixing-length theory (MLT), that describes convective mixing at the radius r close to the 
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boundary. In this model / is a free parameter that is calibrated for each type of connec- 
tive boundary either semi-empirically through observations, or through mult i- dimensional 
hydrodynamic simulations. 



The MESA CBM model is based on the findings in hydrodynamic models (IFreytag. Ludwig. fc Steffen 



19961 ) that the velocity field, and along with it the mixing expressed in terms of a diffusion 
coefficient, decays exponentially in the stable layer adjacent to a convective boundary. Fol- 
lowing these findings, the CBM model extends time-dependent mixing according to the MLT 
diffusion coefficient -Dmlt across the Schwarzschild boundary with the diffusion coefficient 
given by Equation (j2J). The total diffusion coefficient is therefore D = -Dmlt + Adv- I n the 
CBM model adopted here the energy transport in the convectively stable layer is assumed 
to be due to radiation only. Th is CBM model was first introduced in stellar evolution cal- 
culations by lHerwig et al.l ( 119971 ). This, or very similar models, have been applied to several 
related situations in stellar evolution. The most relevant, because similar, case is CBM at 
the bottom of the He-shell flash (or pulse-driven ) convection zone (PDCZ ) in AGB stars (e.g. 
Herwig et al.lll999t iMiller Bertolami et al.ll2006l ; IWeiss fc Fergusonll2009l ). The consequences 
include larger 12 C and 16 abundances in the i ntershell, in agreement with observations of 
H-deficient post-AGB stars if /pdcz ~ 0.008 (IWerner fc Herwigi 120061 ) . Multidimensional 
hydrodynamic simulation s of He-shell flash convection seem to support this value of /pdcz 
( IHerwig et al.ll2006l . 120071 ). but more sophisticated numerical hydrodynamics work is needed. 



The MESA CBM model is meant to represent a wide range of hydrodynamic insta- 
bil ities that may contribute to mixi ng at the convective boundary. The cases considered 
by IFreytag. Ludwig. fc Steffenl ( 119961 ) featured shallow near-surface convection zones with a 
small ratio of the stability in the unstable and stable zones. These convection zones display 
the classical overshoot picture in which coherent convective systems cross the convective 
boundary and then turn around due to buoyancy effects. The boundaries of shell-flash con- 
vection, such as those in novae or in AGB stars, are much stiffer, and coherent convective 
blobs cannot cross the convective boundary. Instead, shear motion, induced by convective 
flows and internal gravity waves lead to mixing at the convective boundaries in which the 



Kelvin - Helmholtz instability plays an important role ( IHerwig et al.l 120061 ; ICasanova et al. 



201 lbl ). The amount of CBM, as expressed in the free parameter /, depends on the details 
of the specific conditions, including the relative stability of the stable to unstable side of 
the boundary as well as the vigour of the convection. While for shallow surface convection 
/ was found to be in the range 0.25 . . . 1.0, hydrodynamic and semi-empirical studies show 
that / = 0.008 is appropriate for the bottom of the He-shell flash convection zone. 

For the nova simulations we adopt /nova = 0.004 at the bottom of the TNR convective 
zone, that eventually includes most of the accreted envelope. This number is of the same 
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order of magnitude (within a factor of two) compared to CBM efficiencies that have been 
found to reproduce observables related to the He-shell flash convection zone in AGB stars. 
As a result, 12 C and 16 are mixed from the WD below, and their combined mass fraction 
in the convective zone reaches Zqq ~ 0.29. Such a CO abundance is similar to the average 
mass fraction of CNO elements observed in the ejecta of CO novae. With this increase in 
the total CO abundance, the simulation produces a fast CO nova. The surface (bolometric) 
luminosity increases by six orders of magnitude on a timescale of 50 seconds (the solid red 
curve in the upper-right panel in Fig.[6]). The radius increases on a longer timescale, of the 
order of 10 4 seconds (the lower- left panel). This corresponds to a surface expansion velocity 
of nearly 300km/s. The velocity exceeds the speed of sound only in the outer layers of the 
expanding envelope, which has a negligible relative mass (the lower-right panel). Given the 
very short evolution timescale of this nova model, its post-TNR mass-loss cannot anymore 
be caused by the super-Eddington luminosity alone because the associated mass-loss rate 
([I]) is too slow. Instead, the envelope would rather quickly fill the WD's Roche lobe, after 
which it would probably be expelled from the binary sys tem as a result o f its (common- 



envelope) interaction with the secondary component (e.g., iLivio et al.lll990l ). The details of 
this process are not yet understood, and we do therefore not attempt multicycle simulations 
of nova models with CBM. 

As an additional test, we compare the maximum hydrogen-burning luminosities Lh 
achieved in the basic 1.2 M CO nova models with simple estimates based on the approxi- 
mation of the nuclear energy generation rate limited by the mass fraction of CNO elements 
Zcno i n the convective envelope of a CO nova during its TNR. When limited by the /3 decays 
(e.g. "Hot" CNO cycle), 

c r 1 nl3 ( ^CNo\ _1 _i 

Uax ~ 5.6 x 10 id I -^j- I erg • g • s \ 



(IGlasner. Livne. &: Truranl |2007| ) . and given that L H ~ £ m ax^conv, where M conv is the mass 



of the convective envelope, the hydrogen-burning luminosity can be estimated as: 

i -^H . ZcNO -Mconv /„\ 

log 10 - « 8.48 + log 10 — + log 10 IF ^-. (3) 



Without any CBM the CNO abundance is Z CNO = Z Q = 0.019 and M conv w M acc = 
1.9 x 1(T 5 M (data from Table □ for M WD = 1.2 M , T WD = 3.0 x 10 7 K, and M = 
lCr 10 M o /yr), in which case the equation fl3]) estimates log^L-^/ L Q ) ~ 9.04, while our 
numerical simulations give 8.74. In models with CBM the mass of the convective envelope 
M conv ps 4.2 x 10~ 5 M exceeds the accreted mass by the amount of material dredged up from 
the WD core. In this case, Z CNO ~ ^co ~ 0.29, which results in log 10 (L H /£o) = 10.56. The 
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numerical simulations give log lo (LH/£ ) = 10.67. The maximum H-burning temperature 
reached during the outburst of our mixed CO nova is 2.32 x 10 s K which is 7 x 10 6 K higher 
than T max in the corresponding unmixed model. 



5.2. Mixing Caused by 3 He Burning 



Shen fc Bildstenl (120091 ) have quantified the role of 3 He in the onset of a nova. They 
have shown that if the mass fraction of 3 He in the H-rich material accreted onto a WD 
is higher than X( 3 He) = 2 x 10~ 3 then convection in the nova envelope is triggered by 
the 3 He( 3 He,2p) 4 He reaction, rather than by 12 C(p,7) 13 N. This alters the amount of mass 
that is accreted prior to a nova outburst and should therefore be taken into account when 
comparing the observed and theor etically predicted nova rates. As a likely place for novae 
with the 3 He-triggered convection, IShen fc Bildstenl (120091 ) consider binaries in which a low- 
mass MS component has undergone such significant mass-loss that it now exposes its formerly 
deep layers where 3 He was produced in a large amount as a result of incomplete pp-chain 
reactions (the so-called " 3 He bump"). 

In our simulations, we have found a variation of this 3 He-triggered convection scenario, 
i.e. that the nova can generate the 3 He in situ. The new scenario is based on the results 
obtained by lTownsley &: Bildstenl (120041 ). who have demonstrated that WDs accreting with 
rates M < 1O _1O M / yr should maintain their central temperatures at the level of Twd < 
10 7 K. We have presented such a CO WD model with M WD = 1.2 M , T WD = 3.3 x 10 6 K 
and M = 10~ n M Q / yr in SectionHl In this model a large amount of 3 He, X( 3 He) « 5 x 10~ 3 , 
is produced at the base of the accreted envelope (the dot-dashed blue curve in the middle 
panel in Fig.JTj). The "sloped 3 He enhancement" is formed as a result of incomplete pp- 
chain reactions, like the 3 H e bump in low-m a ss MS stars. Eventually, 3 He ignition triggers 
convection, as predicted by lShen &: Bildstenl (120091 ). When we include CBM in this model, 
as we did for the standard mixing model using / nova = 0.004 in the equation (J2j) , the 3 He- 
driven convection penetrates into the WD's outer layers and dredges up large amounts of C 
and O into the convective envelope (the middle panel in Fig.[8|), allowing subsequently for 
a fast nova. Note that this interface mixing occurs before the TNR, when the maximum 
temperature is still lower than 5 x 10 7 K. It is not until the 3 He abundance in the convective 
zone decreases below its solar value that the major TNR driven by H-burning in the CNO 
cycle will ensue. 

A comparison of the time intervals between successive nova outbursts (the accretion 
times), t acc = M acc /M, that can be estimated using the corresponding numbers from Tabled! 
for the 1.2 M models with Twd = 20 x 10 6 K and Twd = 3.3 x 10 6 K, and with the accretion 
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rates 10 -10 M /yr and 10 -11 M /yr, respectively, shows that for one event involving the cold 
and slowly accreting WD there should be nearly 59 events occurring on the hot and faster 
accreting WD, provided that the both types of cataclysmic variables are already present 
in equal amounts. Given that the last assumption is not true, because the cooling time 
for the second model is much longer than that for the first one, the relative observational 
frequency of novae with the 3 He-triggered convection should actually be very low. We plan 
to study nova models with the 3 He-triggered convection in more detail, even though as a 
purely theoretical case, in our future work. 



6. Conclusion 



We have presented a grid of nova simulations with the one-dimensional stellar evolution 
code MESA, and provided some tests and comparison with works by others, to demonstrate 
that MESA can generate state-of-the-art nova simulations. In addition, we have investigated 
the effect of convective boundary mixing at the bottom of the TNR convection zone. Inter- 
estingly, the CBM efficiency that reproduces observed CO enhancements is of the same order 
of magnitude (within a factor two) compared to CBM efficiencies that have been found to 
reproduce observables related to the He-shell flash convection zone in AGB stars. We do not 
consider this a coincidence, but rather it is likely that in both cases very similar physics of 
CBM is at play, namely shear motion of convective flows and internal gravity waves leading 
to mixing at the convective boundary, in which the Kelvin-Helmholtz instability plays an 
important role. It will be exciting to study both phenomena side by side in the future. 
For example, the one-dimensional simulations perf ormed here are relat ed to the results of 
multi-dimensional simulations of CBM reported by iHerwig et al.l ( 120071 ). The CBM mixing 
parameter / nova = 0.004 turns out to be sufficient to reproduce both the observed heavy- 
element enrichment of CO novae, as well as amounts of CO-rich material dredged up in 
numerical simul ations of the Kelvin-Helmhol tz and other hydrodynamic instabilities in the 
Casanova etaD bold , boilal lb}) . 



novae case [e.g. 



Further, we have studied slow accretion of solar-composition material onto a CO WD 
wi th the central temper a ture T wo = 3.3 x 10 6 K, exploring the original scenario proposed 
by iTownsley fc Bildstenl (120041 ) in which nova outbursts can actually occur under such con- 
ditions. We have found that incomplete pp-chain reactions lead to the formation of a sloped 
3 He enhancement at the base of the accreted e nvelope in this cas e , with the maximum mass 
fraction X( 3 He) ~ 5 x 10~ 3 . As predicted by IShen fc Bildstenl (120091 ) for this high abun- 
dance, 3 He ignites before the major TNR, and this triggers the development of a convective 
zone adjacent to the WD's surface. When complemented with CBM mixing, sufficiently large 
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amounts of C and are mixed into the envelope to produce a fast CO nova. These new results 
may s uggest a more probable scenario, as compared to the one proposed by lShen &: Bildsten 
( 120091 ). for the 3 He-triggered novae. Clearly this aspect of the nova evolution deserves further 
investigation, in spite of the fact that the relative observational frequency of novae with the 
3 He-triggered convection is expected to be very low. 
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Fig. 1. — Generation of CO WD models with masses 0.85 M and 1M using MESA and 
"stellar engineering" tricks (see text). The computations have used the 6M and 8.5 M Q 
pre-MS stars as the initial models (the pre-MS evolution is not shown). The detailed stellar 
evolution computations (the solid blue curves) are followed by those in which the maximum 
opacity is reduced to a small value (the red curves) and the total mass of the star is relaxed 
to Mwd (the solid black and dashed green curves). After that, the maximum opacity is 
restored (the small red circles) and the WD model is given time to relax (the dot-dashed 
cyan curves). Finally, the WD model is cooled off to a desired temperature (the dotted 
magenta curves). 
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Fig. 2.— A snapshot of the evolution of the 1.15 M CO WD with T WD = 3 x 10 7 K, 
accreting the solar-composition material with the rate 1O _1O M / yr. Left panel: a track in 
the Hertzsprung-Russell diagram for the whole nova multicycle evolution (accretion, TNR, 
mass loss). The star symbol shows the model internal structure of which is displayed on the 
right. The three right panels depict internal profiles of various stellar structure parameters 
as functions of the quantity x = log 10 (l — q), where q = M r j (M WD + M acc ) is a relative mass 
coordinate. Here, M r is the mass inside a sphere of the radius r, while M acc is the accreted 
mass. For Mwd ~ M , a value x of the abscissa approximates to the decimal logarithm of 
solar masses located to the right of this coordinate. Upper-right: temperature (solid red) and 
density (dashed blue); middle-right: mass fractions of some isotopes; bottom right: adiabatic 
V a d, radiative V ra d, and actual V temperature gradients (logarithmic and with respect to 
pressure, dot-dashed blue, dashed green, and solid red lines). In the convectively unstable 
region, V rad > V > V ad . 
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Fig. 3. — A comparison of three nova models that started their evolution with the CO 
WDs having the same mass (Mwd — 1.15 M Q ) and with the same accretion rate (M = 
10~ 10 M o / yr). The differences in their evolution tracks (with dashed black lines of constant 
radius, upper-left panel), luminosity curves (upper-right panel), maximum temperatures of 
H burning (lower- left panel), and accreted envelope masses (lower- right panel) are caused 
by their different WD initial central temperatures (T^d = 3 x 10 7 K for the dashed blue 
and solid green curves, and T WD = 1.5 x 10 7 K for the dotted red curve) and the chemical 
compositions of accreted material (the solar composition for the solid green and dotted red 
curves, and a mixture of 70% solar and 30% CO WD compositions for the dashed blue 
curve) . 
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Fig. 4. — A comparison of the hydrogen burning luminosity (Q and Lh, black curves), surface 
luminosity (L, red curves), and maximum temperature (green curves) evolution profiles from 
Ami Glasner's 1.147 M (the upper panel) and our 1.15 M & (the lower panel) CO nova models 
accreting solar-composition material with the same rate, M = 10~ 10 M Q / yr. The difference 
in the rise of the surface luminosity is probably explained by different initial conditions and 
our better modeling of the evolution preceding the TNR. 
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Fig. 5. — Solar-scaled mass- averaged abundances in the expanding envelope of a last model 
from our simulations of a CO nova with Mwd = 1.15 M©, Twd = 15 x 10 6 K, and M = 
2 x lO -lo M /yr. The accreted material was assumed to be a mixture of 50% solar and 
50% WD's core c ompo sitions. Our final abundances agree very well with those presented by 



Jose fc Hernanzl (119981 ) in their Fig. 1 for a CO nova model with similar parameters. 
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Fig. 6. — Results of the 1.2 M CO nova simulation with convective boundary mixing ac- 
cording to Equation [2] (solid red curves). Upper-left panel: HRD; upper-right and lower-left 
panels: very fast increase of the surface luminosity and a bit slower radial expansion; lower- 
right panel: velocity profile (in the units of the sound speed) in one of the last models. For 
comparison, the dashed red curves show the results for the same model but without CBM. 
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Fig. 7. — Snapshot similar to Fig.[2]for a 1.2 M CO WD model with the central temperature 
T\vd = 3.3 x 10 6 K that accretes with the rate M = 10~ n M / yr. Note formation of a 
sloped 3 He enhancement at the base of the accreted envelope (the dot-dashed blue curve in 
the middle panel) during the accretion. The envelope is convectively stable yet, therefore V 
coincides with V ra d everywhere (the lower panel). 
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Fig. 8. — Same as Fig.[7J but at a later time. 3 He ignition triggers a convective zone (the 
region where V ra d > V > V a d in the lower panel). Note the enhanced amount of C and O 
(middle panel) in this model with CBM, sufficient to produce a fast CO nova. 
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Table 1. Characteristics of Our CO Nova Models 
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a This model accretes a mixture of 70% solar and 30% WD's compositions. 
b This simulation includes CBM © with / = 0.004. 



